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ABSTRACT 

We study the properties of a sample of 211 heavily-obscured Active Galactic Nucleus (AGN) can- 
didates in the Extended Chandra Deep Field-South selecting objects with /24/im//i?>1000 and R- 
K>4.5. Of these, 18 were detected in X-rays and found to be obscured AGN with neutral hydrogen 
column densities of ~10 23 cm -2 . In the X-ray undetected sample, the following evidence suggests a 
large fraction of heavily-obscured (Compton Thick) AGN: (i) The stacked X-ray signal of the sample is 
strong, with an observed ratio of soft to hard X-ray counts consistent with a population of ^90% heav- 
ily obscured AGN combined with 10% star-forming galaxies, (ii) The X-ray to mid-IR ratios for these 
sources are significantly larger than that of star-forming galaxies and ~2 orders of magnitude smaller 
than for the general AGN population, suggesting column densities of A^j>5xl0 24 cm" 2 , (iii) The 
Spitzer near- and mid-IR colors of these sources are consistent with those of the X-ray-detected sample 
if the effects of dust self- absorption are considered. Spectral fitting to the rest-frame UV/optical light 
(dominated by the host galaxy) returns stellar masses of ~10 11 M Q and <E(B-V)>=0.5, and reveals 
evidence for a significant young stellar population, indicating that these sources are experiencing con- 
siderable star-formation. This sample of heavily-obscured AGN candidates implies a space density at 
z~2 of ~10 -5 Mpc -3 , finding a strong evolution in the number of Lx>10 44 erg/s sources from 2=1.5 
to 2.5, possibly consistent with a short-lived heavily-obscured phase before an unobscured quasar is 
visible. 

Subject headings: galaxies: active, Seyfert; X-rays: galaxies, diffuse background 



1. INTRODUCTION 

Understanding the processes of galaxy formation and 
evolution is one of the outstanding problems of modern 
astronomy. It is now clear that the growth and feeding of 
the central black hole of an active galactic nucleus (AGN) 
and the formation of its host galaxy must be intimately 
connected, as evidenced by the striking correlations be- 
tween the properties of supermassive black holes and 
the stellar systems in which they reside jM agorria n et al.l 
[1991: [Gephardt et al.ll2000l : iFerrarese et al.ll200lD . Fur- 
thermore, it is now strongly suspected that AGN activity 
is critical in the regulation of star format ion in the host 
galaxy either by removi ng all the gas fe.g. JHookins et all 
120061 : [Menci et al.ll2006h . by heating it (e.g. ICroton et all 
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120061 ) or both (jSchawinski et al.l l2006f K However, an 
observational link between AGN activity and quench- 
ing of star formation h as not yet been established (e.g., 
iSchawinski et alj [20091 ). To understand how the AGN 
and galaxy formation processes are connected requires a 
detailed census of the AGN population as well as detailed 
study of individual objects. 

The visibility and detectability of an AGN is deter- 
mined by the interplay among black hole mass, accretion 
rate, gas distribution, dust geometry and the amount of 
star formation in the host galaxy. Thus, no single selec- 
tion method provides a complete view of how galaxies 
and AGN co-evolye. M an y surveys (e.g.. iBarger et aD 
120011: ICowie et al.ll2002l : see lBrandt fc Hasingerll2005l for 
a review) have focused on X-ray selection because X-ray 
emission is a universal feature of AGN accretion. How- 
ever, X-ray selection has limited sensitivity to heavily 
obscured and /or low accretion rate s ources, even in the 
deepest surveys ([Treister et al.l 12004). In particular, ex- 
amples of the most heavily-obscured AGN, the so-called 
Compton-Thick (CT) AGN, which have r-1 for Comp- 
ton scattering or A^^IO 24 cm" 2 , are mostly found in 
the local Universe. Only a few of these sources have 
been identified at higher re dshifts, although this is sim- 
ply due to selection effects ([Delia Ceca et al1l2008bl) . In 
contrast, early AGN population synthesis models that 
can fit the X- ray background (XRB) spectral shape 
and intensity ([Treister fc Urrvl 120051 : iGilli et all |2007| ) 
require ~20-30% CT AGN in order to explain the ob- 
served peak at ~30 keV in the XRB radiation. Re- 
cently, Trei ster et al.l (|2009a[ ) combined observations of 
local (z<0.05) AGN at very hard X-rays (E>20 keV), 
using the International Gamma-Ray Astrophysics Lab- 
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oratory (INTEGRAL) and Swift satellites, with XRB 
models and concluded that the CT AGN fraction is only 
~10% at these redshifts. They conclude that CT AGN 
represent a similarly small fraction of the total XRB ra- 
diation and total black hole growth. However, because 
only ~2% of the XRB comes from CT AGN at z>2 and 
only a few sources are known at such high redshifts, this 
population still remains basically unconstrained, with up 
to an order of magnit ude uncertainty in number density 
(jTreister et al.ll2009aD . 

Because the energy absorbed at optical to X-ray wave- 
lengths is later re-emitted in the mid-IR, it is expected 
that AGN, in particular the mo st obscured ones, sh ould 
be very bright mid-IR sources (jTreister et alj 120061 and 
references therein) , easily detectable in observations with 
the Spitzer observatory. Sources having mid-IR ex- 
cesses, relative to their rest-frame optical and UV emis- 
sion, have been i dentified as poten ti al CT AGN can- 
didates at z~2 (iDaddleLa lJ l^Og; iFiore et alJ 120081 : 
iGeorgantopoulos et al.l [20081 and others). However, be- 
cause of the strong connection between vigorous star for- 
mation and AGN actiyity in the most luminous infrared 
sources ([Sanders et alJ [i~988) , the relative contribution of 
these pro cesses is still uncer t ain and remains controver- 
sial (e.g- lDonlev et al.ll2008t iPope et al.ll2008h . 

In order to detect CT AGN up to high redshifts, 
deep multiwavelength coverage is critical. The Extended 
Chandra Deep Field South (ECDF-S) is one of the best 
fields to carry out this study. The rich multiwavelength 
data available in the ECDF-S includes four ^250 ksec 
Chandra pointings covering an area of ~0.3 deg 2 , and 
deep Spitzer data in the IRAC bands from the Spitzer 
IRAC/MUSYC Public Legacy in ECDF-S (SIMPLE) 
survey and at 24 fim from the Far-Infrared Deep Ex- 
tragalactic Legacy Survey (FIDEL). The 24 /im depth in 
this field is ^35 /iJy, comparable to the deep Great Ob- 
servatories Origins Deep Survey (GOODS) observations. 

In this paper, we present a study of the properties of 
the CT AGN candidates in the ECDF-S selected from 
their 24 /im to optical flux ratio and rest-frame optical 
to UV colors. In Section 2 we describe the multiwave- 
length data used in this work, while our selection scheme 
is presented in Section 3. The properties of the sources 
individually detected in X-rays are discussed in Section 
4, while the remaining sources are studied in Section 5. 
The near and mid IR colors, X-ray to mid-IR flux ratios 
and space density of our sources are presented in Section 
6, and our conclusions are reported in Section 7. We 
assume a ACDM cosmology with h —0.7, fl m —0.3 and 
fiA=0.7, in a greement with the m ost recent cosmological 
observations (jSpergel et al.ll2007h . 

2. MULTIWAVELENGTH DATA 

2.1. X-ray Data 

The full ECDF-S was covered by Chandra as part 
of a Cycle 5 guest observer program (PI: N. Brandt). 
The total area covered is ~0.3 deg 2 to a depth of ~230 
ksec. Details about these observati ons, images and cata - 
logs were present e d sep arately by ILehmer et al.1 (120051) 
and IVirani et alj ([20061) . The work of ILehmer et all 
( 2005) reports the finding of 762 sources to flux limits 
of 1.1 xl0~ 16 and 6.7xl0" 16 erg cm^s" 1 in the soft 
(0. 5-2 keV) an d hard (2-8 keV) bands. The catalog 
of IVirani et al.l (|2006h includes 651 sources to similar 



flux limits, using a more conservative rejection of pe- 
riods of higher background, thus reducing the number 
of spurious sources but excluding a few real sou r ces as 
well. In this work, we use mainly the IVirani et~aTI (p006T ) 
catalog. However, when performing stacking analysis 
we further exclude sources individually detected in the 
ILehmer et al.1 (|2005l ) catalog and/or in the deeper 1 and 2 
Msec obse rvations available in the centr al c^O.l deg 2 , pre - 
sented by lAlexander et all (|2003l ) and ILuo et all (|2008l ) 
respectively. 

2.2. Mid-IR Data 

The ECDF-S was observed extensively at near and 
mid-IR wavelengths by the Spitzer Space Telescope. 
The central 10'xl6' region was covered using bot h the 
Infrared Array Camera (IRAC; iFazio et al.l I2004T) and 
the Multiband Im aging Photometer for Spitzer (MIPS; 
iRieke et "all l2004f ) . These observations were performed 
as part of the GOODS survey and guaranteed time 
programs. More details about th ese observations were 
presented by iTreister et al.l (|2006h and R. Chary et al. 
(2009, in prep.). 

The extended 30'x30' region was observed by IRAC 
as part of the SIMPLE survey. The flux limits for 
the SIMPLE observations are 0.76, 0.4, 5.8, and 3.6 
/iJy in the 3.6, 4.5 and 5.7, and 8 /im bands, respec- 
tively, or ^3-5 times shallower than the GOODS ob- 
servations. More details abo ut the IRAC cove r age o f 
the ECDF-S were reported by ICardamone et al.l (|2008l ) 
and M. Damen et al. (2009, in prep.). When required, 
the conversion facto rs from flux densit y to Vega mag- 
nitudes provided by IFazio et al.l (|2004f ) were assumed. 
The ECDF-S was also observed at longer wavelengths in 
the FIDEL survey, which obtained data at 24, 70 and 
160 /im. Because of the significant decrease in sensitiv- 
ity at longer wavelengths, mainly due to source confu- 
sion, only the 24 /im data are considered in this work. 
Details about the data reductio n and c atalog creation 
were presented bv ITreister et~aTl ([2009p[ ). The approxi- 
mate flux limit of the 24 /im data used here is ~35/iJy. 
For comparison, the flux limit of t he GO ODS 24 /im 
data is ~12/iJy ([Daddi et al.ll2007h . while IFiore et all 
(2008) studied only sources with /24>40/iJy in the Chan- 
dra Deep Field S outh (CDF-S) proper. The work of 
IFiore et all pfM ) in the COSMOS field used only the 
shallow cycle 2 MIPS data to a flux limit of ~550 /iJy. 

2.3. Optical/Near-IR Data 

The ECDF-S was observed at optical wavelengths us- 
ing both ground-based and space telescopes. Deep im- 
ages (to V^26.5) in the UBVRI were obtained using 
the Wide Field Imager (WFI) on the 2.2m telescope at 
La Silla. These data were made public b y the ESO 
Dee p Public Survey ([Mignano et all I2007D . COMBO- 
17 (IWolfet all 120041) . and Ga rching-Bonn Deep Survey 
(GaBODS: iHildebrandt et al.l l2006h teams. Imaging to 
z'=23.6 AB, was performed using the MOSAIC-II cam- 
era mounted on the B lanco 4m telescope at Cerro Tololo 
([Gawiser et alJl2006bl ) . The Hubble Space Telescope ob- 
served the ECDF-S in the V and z' bands as part of the 
Galaxy Evolution from Morphology and SEDs (GEMS) 
survey (jRix et al.H200l . Deep near-IR coverage was ob- 
tained using the CTIO 4m telescope with the Infrared 
Sideport Imager (ISPI), reaching a magnitude limit of 



Heavily-Obscured AGN in Star-Forming Galaxies 



3 



-22 (AB) in the JHK S bands (|Tavlor et all 120091) , 
as part o f the Multiwayelength Survey by Yale-C hile 
(MUSYC; iGawiser et all[2006al iTreister et aLll200l . 

In order to obtain high-quality photometric redshifts, 
the ECDF-S was imaged in 18 medium-band optical fil- 
ters using the wide-field Suprime Camera on the Subaru 
telescope. These filters range from 4270 to 8560A and 
are optimized for phot ometric redshift determinations 
(Havashin o et alj [2000). The data reach an average 
depth of ~26 magnitudes (5er, AB), extending redshift 
determinations to sources much fainter than most spec- 
troscopic redshift surveys. More details about these ob- 
servations, data reduction and catalogs will be presented 
by C. Cardamone et al. (2009, in prep.) and in Sec- 
tion [231 



2.4. Spectroscopic Redshifts 

Rich optical spectroscopic data exist in the ECDF-S. 
We took advantage of the ESO/GOODS-CDFS spec- 
troscopy master catalogue 11 which compiles spectro- 
scopic data from 11 different sources, mostly using 8-10 
meter class telescopes. In particular, many of the spec- 
troscopic redshifts for t he X-ray sourc e s in t he CDF-S 
proper were reported by ISzokolv et al.1 (|2004f ) using the 
VLT telescopes with the FORS1/FORS2 spectrographs. 
In addition, spectroscopic observations for 339 X-ray 
sources in the ECDF-S field usin g the VLT/VIMOS an d 
Magellan/IMACS are reported bv lTreister et all (|2009bf ). 

Given the very faint optical magnitudes of the ECDF- 
S X-ray sources, obtaining redshifts and identifications is 
a very hard, time-demanding task. In particular, optical 
spectroscopy is often out of reach, even for state-of-the- 
art lOm-class telescopes. However, near-IR spectroscopy 
offers a higher chance of success, not only because these 
sources are significantly brighter at these wavelengths, 
but more importantly because, based on the experi- 
ence with similar sources in other fields, they are found 
at redshifts z=l-3, where the typical optical emission 
lines are shifted to the near-IR. For the brighter sources 
we used the Spectrograph for INtegr al Field Observa- 
tions in the Near Infrare d (SINFONI; Eiscnhauc r et al.l 
120031: iBonnet et all |2004|) at the VLT, and the Multi- 
Object InfraRed Camera and Spectrogr aph (MOIRCS; 
llchikawa et~aTI 120061 ; ISuzuki et all l2008f) at the Subaru 
8m telescope. For the remaining sources in our sample 
we have to rely on photometric redshifts, which are rela- 
tively accurate because of the deep Subaru medium-band 
photometry described below. 

2.4.1. SINFONI Data 

A total of 16 hours were granted to observe four 
X-ray sources in the ECDF-S with V-K>7 using the 
VLT/SINFONI Integral Field Unit (IFU) as part of pro- 
gram 080.A-0612. The observations were carried out in 
service mode by the Paranal staff in the October 2007- 
January 2008 period. We attempted to observe one of 
the sources (XID 647) with the Laser Guide Star Fa- 
cility; however, due to technical problems, these obser- 
vations were cancelled. Of the remaining three sources, 
two have IR-red excesses, XIDs 277 and 580, and one, 
XID 57, was not detected at 24 fim despite a very red 



optical/near-IR color (R-K=6.56), and hence is not in- 
cluded in the sample discussed in this paper. 

For these sources we used SINFONI in the 0.25" spatial 
resolution mode, which provides a field of view of 8"x8". 
The H+K grism, which provides a wavelength coverage 
of 1.45-2.45 /im and a resolution of R~1500, was used. 
The total exposure time was 3 hours for XID 580 and 
2.7 hours for XID 277. In each case, a square dithering 
pattern with 4" offsets was performed so that the source 
was always in the field of view. This pattern was repeated 
four times. 

Data were reduced using the official ESO SINFONI 
pipeline vl.9.4 12 . The SINFONI data reduction cook- 
book 13 , provided by the Paranal Science Operations 
group was followed. Briefly, raw data were first cleaned 
to remove bad lines created by the presence of hot pix- 
els in the bias region. Then data are bias-removed and 
flat-fielded. Wavelength calibration was performed using 
arc lamps. Distortions were corrected using the provided 
wavelength maps. Finally, a data-cube was reconstructed 
using the sinfo_rec_jitter routine. Flux calibration and re- 
moval of telluric features was performed using a standard 
star, typically a solar analog, observed within 2 hours of 
each science observation. Cube visualization and extrac- 
tion of a 2-dimensional pseudo-slit was done using the 
QFitsView tool 14 . Extraction of the 1-dimensional spec- 
trum, smoothing and flux calibration were performed us- 
ing standard IRAF tools. 

2.4.2. MOIRCS Data 

Source XID 322 was observed using the MOIRCS near- 
IR multi-object spectrograph mounted at the Subaru 
telescope on the night of December 18, 2008. MOIRCS 
was used in its multi-object mask mode, in order to use 
a sample of nearby bright stars to do the source acqui- 
sition. The HK500 grism, which provides a resolution 
of ^700 in the 1.3-2.5 /im wavelength range, was used. 
Ten-minute individual integrations were performed, after 
which the source was dithered by 3" along the slit fol- 
lowing a standard ABBA routine in order to allow for an 
accurate sky subtraction. The total exposure time was 
3 hours. Data reduction was done using the IDL-based 
MOIRCSMOSRED package created by Youichi Ohyama. 
Basic steps include extraction of the slit spectrum, sky- 
subtraction, flat fielding, flexure correction, co-adding of 
individual AB spectra to finally combine the resulting 
sky-subtracted spectra. Wavelength calibration was per- 
formed using OH sky-lines as reference and flux calibra- 
tion was based on observations of HD20758, a G2V star, 
taken during the same night. 

2.5. Photometric Redshifts 

In order to obtain highly accurate photometric red- 
shifts for all of the optically detected sources in the re- 
gion of the ECDF-S, we obtained 18 medium band im- 
ages of the field. Our deep medium-band Subaru imag- 
ing reaches Rab ~ 26, 2 magnitudes deeper than previ- 
ous CDFS medium band imaging (COMBO-17 reaches 
Rab ~ 24 with 12 me dium band filters and 10% accu- 
racy; IWolf et aT1l2004[ ). and provides accurate redshifts 



The SINFONI pipeline can 



be 



found 



http: / / www.eso.org/sci / data-processing/ software/pipelines / sinfoni / sinfo-pipe-r< 
1J http:// www.eso.org/ sci / facilities / paranal /instruments / sinfoni / doc/ 



11 Available at http:/ /www. eso.org/science/goods/spectroscopy/CDFS_rJfas4imili40 c at http:/ /www. mpe.mpg.de/~ott/QFits View/ 
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(~ 1% accuracy) and detailed spectral energy distribu- 
tions for > 90% of the detected X-ray sources. The Sub- 
aru images reach nearly all of the X-ray selected AGN. 
The greater number of filters and their even spacing 
give an effective spectral resolution of R — A/AA ~ 23 
(|Taniguchil 120041 ). The catalog was created using th e 
deep BVR image for detection (|Gawiser et al.l l2006al) . 
and contains roughly 60,000 sources with detailed spec- 
tral energy distribution. The catalog includes color- 
matched photometry from all 18 medium bands, in ad- 
dition fluxes from the optical an d near-IR ground based 
imaging (see iTavlor et al.l 120091 for more details on the 
broad band data) . Details of the medium-band data re- 
duction and characteristics are provided in Cardamone 
et al. (2009, in prep). 

In order to obtain accurate photometric redshifts, we 
used EAZY (Easy and Accurate Zphot from Yale), a pro- 
gram op timized to prov ide high quality redshifts over 
0<z<4 (jBrammer et al.l l2008h . EAZY is a new user 
friendly, full featured redshift fitting c ode, including a 
user- friedly interface based on Hyp er Z dBolzonella et al.l 
120001) and the use of priors (e.g., iBemtej |2000|) . The 
template set and the magnitude priors are based on 
semi-analytical models which are complete to very faint 
magnitudes rather than highly biased spectroscopic sam- 
ples, and so are particularly useful for samples of ob- 
jects such as our optically- faint mid-IR-selected sources, 
where complete spectroscopic calibration samples are not 
available. Further, EAZY introduces a template error 
function to account for wavelength-dependent template 
mismatch rather than relying on minizing the scatter be- 
tween photometric and spectroscopic redshifts of the sub- 
sample of sources bright enough to have spectroscopic 
redshifts. EAZY allows for linear combinations of the 
templates sets. We make use of this feature including 
an additional broad-line AGN template when fitting the 
X-ray sources. EAZY assigns each object a redshift by 
marginalizing over the full posterior probability distribu- 
tion rather than maximizing the likelihood. 

We compute the photometric redshifts for the sources 
detected in the combined BVR image using the default 
EAZY parameters, adopting the default template error 
function and R-band photometric prior, and including 
the optical and n ear-IR ground base d coverage from the 
MUSYC survey (j Taylor et al.l I2009D in addition to our 
medium band data. We find a median Az/(1 + z) of 0.01 
for all sources with spectroscopic redshifts out to z~5. 
Limiting ourselves to the subsample of X-ray sources, we 
do only slightly worse, Az/(l+.z)=0.011, but the fraction 
of outliers doubles to 20%. Overall, for Az/(l+z), this is 
a factor of t hree better than is curr ently done using broad 
band data (jBrammer eLalJ |2008|) : We find good qual- 
ity photometric redshifts (q z < 1; iBrammer et al.l [2008T) 
down to R-band magnitude of AB=26 for 60% of the 
sample. For many of the remaining sources poor fits are 
obtained because the photometry is too uncertain, due 
to faintness of the sources, or intrinsic variability in the 
source over the time period of which the photometry was 
taken ([Salvato et al.l l2009h . There are also cases where 
the intrinsic SED may not be matched by those provided 
in the template set or degeneracies in color-z space can 
result in multiple peaks in the redshift-probability distri- 
bution. If no restriction in the quality parameter, q z , is 
used the statistical quality of the photometric redshifts 



degrades only slightly, to Az/(l+z)=0.016 with an out- 
lier fraction of 27%. Overall, with medium band pho- 
tometry we have achieved highly accurate photometric 
redshifts for the majority of the sources in the ECDF-S 
down to very faint magnitudes. 

3. OPTICAL/MID-IR SELECTION OF CT AGN 
CANDIDATES 

X-ray observations have been very efficient in finding 
unobscured and moderately obscur ed AGN up to high 
redsh ifts and low luminosities (e.g., iBrandt fc Hasingerl 
2005). However, highly obscured, particularly Compton- 
Thick AGN, have been significantly excluded even from 
the deep Chandra and XMM-Newton surveys. Because 
most of the radiation absorbed at X-ray and UV wave- 
lengths is then re-cmittcd in the mid- and far-IR, re- 
cent studies at these energies, mostly taking advan- 
tage of Spitzer observations, have been very successful 
in finding heavily-obscured AGN ca ndidates missed by 
X-ray selection up to h igh redshifts (IDaddi et alj 120071: 
Alexander et~afl l2008t iDonlev et all 120081 : iFiore et all 
20081 ). 

In particular, IFiore et al.l (|2008l ) presented a selec- 
tion method based on the 24 fim to R band flux ratio 
and R-K color; specifically, /24///i>1000 and R-K>4.5 
(Vega) . According to simulations based on stacking of 
the X-ray signal, they estimate the fraction of heavily- 
obscured AGN in this sample to b e greater than 80% 
(jFiore et al. 1 1 20081 120091 ) . Similarly, IDaddi et all ([20071 ) 
found a large number of CT AGN candidates by select- 
ing sources which show a significant excess in the star 
formation rate measured from IR light, co mpared to the 
value derived in the ultraviolet. However. iMurphv et al.l 
(120091) showed, using mid-IR spectroscopy and adding 
far- IR photomet r ic dat a for the sources selected using 
the IDaddi et all (|2007T ) method, that in only -30% of 
the cases the observed mid-IR excess was due to the pres- 
ence of an AGN, while in the majority of the sources the 
presence of strong aromatic PAH features explained the 
discrepancy in the derived star formation rates. 

We applied the "infrared-red" selection criteria of 
IFiore et al] (|2008l ) to the MlPS-selected sources in the 
ECDF-S to select CT AGN candidates, as shown in 
Fig. □ Of the 7201 24/mi sources detected in the FI- 
DEL observations to a flux limit of —35 /LtJy, 211 (—3%) 
satisfy the ./V/aMOOO an d R-K>4.5 cuts. O f the 651 
X-ray detected sources in the lVirani et al] (120061 ) catalog, 
18 are found in this region, —2.8%, a similar fraction as 
in the general population. The fraction of infrarcd-rcd 
sources with a direct X-ray detection i s —9%, smaller 
than the 16% reported by IFiore et al.l (|2008f) for the 
CDF-S observations, and significantly smaller than the 
fractio n of direct X-ray detections in the COSMOS field 
(40%: IFiore etal] 120091 . As can be seen in Fig. [TJ the 
typical value of /24//.R for the ECDF-S sources is —50- 
100, while <i?-A'>— 3.5, so the IR excesses are extreme, 
making them good candidates to be heavily obscured 
AGN. 

4. X-RAY DETECTED AGN 

A total of 18 sources in the IR-red excess region were 
significantly detec ted individua l ly in X -rays and included 
in the catalog of IVirani et alj p006). In addition, six 
more sources were detected in the 2-Msec C handra ob- 
servations covering the central CDFS region (|Luo et al.l 
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Fig. 1. — Ratio of 24 /im to R-band flux density as a function 
of R-K color. Contours show the location of all the MlPS-selected 
sources. Red circles show the location of X-ray sources with 
HR (define as (H-S) / (H+S) where S and H are the background- 
subtracted counts in the soft and hard bands respectively) greater 
than -0.3, while blue circles have HR<-0.3. Small black circles 
show the 193 sources in the IR-red excess region that were not de- 
tected individually in X-rays. While there is no clear trend of HR 
with optical (R-K) color, most of the X-ray hard sources also have 
large R-K colors. 



2008) ; two of the l atter sources were included in the 
Lehmer et al.l (120051) catalog, but were not reported by 



Virani et alj ( 20061). Four of these six sources were not in- 
cluded in the lAlexander et alj (j2063) catalog, indicating 
that they were very faint in X-rays and only significantly 
detected with the a dditional 1 Msec o f data. 

For comparison, iFiore et al.l (|2008f ) report that the 
same number of IR-red excess sources, 18, were directly 
detected in X-rays and included in the 1 Msec Chandra 
catalog, with four more marginally detected. This can be 
expected by the combination of deeper X-ray data (by a 
factor of 4) and smaller area (by a factor of ~7). Also, 
an important factor is the difference in depth of the K 
band image used in their analysis. While the MUSYC 
ECDF-S images reach a magnitude limit of if ~20.2, the 
GOOD S-MUSIC catalogs r each almost two magnitudes 
deeper ()Grazian et~aT][200l . 

After an extensive archival search, we found opti- 
cal spectroscopy for two sources in our X-ray-detected 
sample, namely XIDs 480 and 500. Source XID480 
is an intriguing so urce. Its redshift, z= 1.603, was se- 
curely measured bv lSzokolv et al.l (|2004l ) based on deep 
VLT/FORS1-FORS2 observations. These authors clas- 
sified XID 480 as a QSO-1, based on the soft X-ray spec- 
trum. However, only narrow emission lines are present 
in the optical spectr um. From X-ray spectral fitting 
(jTreister et al.ll2009bf) . we found that N H =2x 10 22 cm~ 2 , 
consistent with the observed hardness ratio. Hence, this 
source is at the unobscured/obscured AGN boundary, 
which explains the discrepancy between the observed 
soft X-ray spectrum and the red optical/near-IR colors. 
Source XID 500 was observed by VLT/V IMOS as part 
of our MUSYC identification program (|Treister et al.l 
2009b). The optical spectrum of this source is shown 
in Fig. [2] The spectroscopic redshift for this source is 
3.343, based on a strong narrow emission line at 5402 
A identified as Lya. In support of this interpretation, a 
strong decrement in the continuum is observed blueward 
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Fig. 2. — Opt ical spectrum of X- ray source 500 obtained using 
VLT/VIMOS UTreister et al.H2009bfl . Wavelength is shown in the 
observed frame. The redshift for this source is 2=3.343, identify- 
ing the strong emission line at ~5400 A as Lyo, which is sup- 
ported by the significant continuum decrement bluewards of this 
line. The photometric redshift for this source is 1.75±0.1 (99% 
confidence level), but with a very bad x 2 - The closest identifica- 
tion matching the phometric redshift and the observed emission 
line is [CIII]1909A at z=1.83; however, if this is the case, CIV at 
observed-frame 4383 A should also be visible. Hence, most likely 
the photometric redshift is wrong in this case. 



of this line, while the continuum is clearly detected on 
the red side. 

The near-IR spectra of sources 277 and 580 can be 
seen in Fig. [3l The redshift of XID 277 is secured by the 
presence of a strong emission line, which is marginally 
resolved at the SINFONI resolution and consistent with 
the wavelengths of Ha and [Nil] at a redshift of 1.286. 
In addition, a weak [SII] line is present in the data. For 
XID 580, no obvious emission line is detected. However, 
on further inspection, two weak emission lines are found 
at wavelengths of 1.866 and 2.08 microns. These lines 
could be identified as 7Ja+[NII] and Hel respectively at 
a redshift of 1.845. The weak detection of these lines can 
be explained by their location in a spectral region with 
very low atmospheric transmission. The main feature in 
the near-IR spectrum of XID 322 is a strong emission 
line 1.722 microns. At the resolution of our setup, this 
line is barely resolved into two peaks, the largest one at 
1.7179 /im and a smaller one at 1.7257 microns. This is 
consistent with Ha and [Nil] at a redshift of 1.621. No 
other feature is visible in this spectrum. 

The redshift distribution for the IR-red excess sources 
detected in X-rays, including both photometric and spec- 
troscopic measurements, is shown in Fig. 01 The aver- 
age redshift for our sample is 2.37 (m edian=1.85) , larger 
than the value of 1.55±0.53 found bv lFiore et all (|2009l ) 
for similar sources in the COSMOS field. This is most 
likely due to the much brigh t er 24 j jltcl flux limit of their 
sample. Indeed, IFiore et alJ (|2008l ) reported an average 
redsh ift of 2.1 in the deeper CDF S observations. Simi- 
larly, [GiorgMto20uloi^Fal] (|2008[ ) estimated an average 
redshift of ^2 for similar sources in the Chandra Deep 
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Fig. 3. — Near-IR spectra for three X-ray-detected IR-red ex- 
cess sources in the ECDF-S. Wavelength is shown in the observed 
frame. Shaded wavelengths show spectral regions affected by low 
atmospheric transmission and OH sky lines (Roussclot ct al. 2000). 
These spectra were obtained from VLT/SINFONI (XID-277 and 
XID-580) and Subaru/MOIRCS (XID-322). 




Fig. 4. — Left panel: Redshift distribution for the 24 /mi-selected 
sources in the ECDF-S. The solid histogram shows the distribution 
for the sources not detected in X-rays, while the dashed hatched 
histogram considers only the X-ray detected sources. A KS test 
shows that the hypothesis that these two distributions were drawn 
from the same parent distribution is refuted at the 84% confidence 
level. The dotted histogram shows the slightly lower redshift distri- 
bution (divided by 1,000) for all the sources with a 24 (im detection 
and a measured photometric redshift in the ECDF-S. Right panel: 
distribution of R-band magnitude for the IR-red excess sources de- 
tected in X-rays (solid histogram), those not detected in X-rays 
(dashed histogram) and all sources with a measured photometric 
redshift (hatched histogram). Not surprisingly, spectroscopic and 
photometric redshifts are available preferentially for the optically- 
brightest sources. In addition, X-ray detected sources are brighter 
in the optical than the non- X-ray detected ones. 



Field North. For our sample, the minimum redshift is 
1.286 (XID 277) and the maximum is 4.65 (XID 303), 
although the latter is only a photometric redshift. 

4.1. X-ray Properties 

Automated X-ray spectral fitting was performed for 
the sources in the ECDF-S using the Yaxx 15 software. 
An absorbed power-law was assumed for all sources. 
For the two brightest sources, XID 480 and 284, with 

15 Available at http://cxc.harvard.edu/contrib/yaxx/ 
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Fig. 5. — Neutral hydrogen column density (Nh) as a function 
of hardness ratio for the 15 X-ray sources with either spectroscopic 
or photometric redshifts. Njj was calculated from spectral fitting 
for the X-ray brightest sources and from the HR for the remaining 
ones. For the two brightest sources, XID 480 and 284 we measured 
r=2.5±0.6 and 1.7±0.7 respectively. For the remaining sources, 
we assumed an intrinsic power-law spectrum with T=1.9 at the 
observed redshift of the source and correcting for galactic absorp- 
tion. In general, these are very obscured sources in X-rays, with 
N[{>10 22 cm -2 , and two have Compton-thick absorption levels 
(XID 303 and 585). 



more than 200 background-subtracted counts, the spec- 
tral slope and amount of absorption were fitted simul- 
taneously. In both cases, the observed sp ectral slope is 
consi stent with the fiducial T=1.9 value (|Nandra et al.l 
Il997h . For fainter sources, the spectral slope was fixed 
to r=1.9 and only the amount of absorption was fit- 
ted. For sources in which the number of X-ray counts is 
even smaller (<80 counts), it is not possible to perform 
spectral fitting. In those cases, we use the hardness ra- 
tio (HR), defined as the ratio between the difference in 
counts in the hard and soft counts and their sum; the 
amount of absorption can be estimated from the HR by 
comparing the observed value with the expected counts 
for an intrinsic obscured power-law with T=1.9 at the 
redshift of the source and folding in the response of Chan- 
dra/ACIS. Because the same intrinsic spectral slope is 
assumed by the HR and spectral fitting methods, consis- 
tent results are obtained. The observed HR values and 
the corresponding calculated Nh for our X-ray detected 
sample are shown in Fig. [5l 

From our derived Nh values we found two sources 
(13%) with iV H <10 22 cm" 2 , i.e., unobscured AGN, 11 
(73%) with 10 22 <iV ff <10 24 cm" 2 (moderately obscured 
sources), and two Compton-thick AGN. The ratio of ob- 
scured to unobscured AGN in this sample, 13:2, is sig- 
nificantly higher than the average value o f ^3-4:1 found 
in de ep surveys with no color cut (e.g., iTreister et al.1 
I2004D . Furthermore, all the sources in this sample have 
Lx>10 43 erg s _1 . Due to the strong depen dence of the 
fraction of obscured AGN on luminosi ty (e.g .. lUeda et all 
l2003tlSteffen et al.ll2003HBarger et al.ll2005fi . the fraction 
of obscured to unobscured sources is expected to be sig- 
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nifican tly lower (^1:2 at this luminosity; iTreister et al.l 
2009b), making the excess in the relative number of ob- 
scured AGN more significant. This confirms that the IR- 
red excess selection preferentially find obscured sources, 
although only a small fraction of this sample reaches 
Compton-thick absorption levels. 

5. X-RAY UNDETECTED CT AGN CANDIDATES 

As was mentioned in Section [3l a total of 211 24-/im 
sources in the ECDF-S were found in the IR-red ex- 
cess region, 18 of them individ ually detected in X -rays 
and included in the catalog of IVirani et al.l (|2006h . Of 
the re maining sources , eight were detected in X-rays ei- 
ther bvlLehmer et alj (|2005[ ) or in the deeper catalogs of 
I Alexander et all (|2003l )~a"nd ILuo etldl (|2008l ). Here we 
focus on the 185 X-ra y undetected obscured-AGN can- 
didates. According to iFiore et al.l ([2008) the vast ma- 
jority of these sources should be heavily obscured, even 
Compton-thick, AGN, in which the combination of pho- 
toelectric absorption and Compton scattering explains 
the lack of X-ray detection. In order to test this hypoth- 
esis, we study the nature of these sources at many wave- 
lengths through X-ray stacking, optical/mid-IR spectral 
analysis and IRAC colors. 

5.1. Photometric Redshift Distribution 

Contrary to the situation for the X-ray detected 
sources, where four sources had a measured spectroscopic 
redshift, none of the 185 X-ray undetected sources have a 
spectroscopy redshift. Hence, redshifts for these sources 
will be based on photometric measurements only. Be- 
cause of the faintness of the optical/near-IR counter- 
parts, accurate photometric redshifts are available for 
only 90 (49%) sources. The remaining 95 sources were 
not detected in the deep BVR image and hence it was not 
possible to compute photometric redshifts. The redshift 
distribution of these sources is shown in Fig. |4] 

The vast majority of the mid-IR selected sources have 
photometric redshifts greater than one, with a median of 
1.9, a significantly smaller value than for the X-ray de- 
tected sources. Performing a Kolmogorov-Smirnov (K-S) 
test to these two distributions we found that the hypoth- 
esis that they were drawn from the same parent distri- 
bution is marginally rejected at a 84% confidence level. 
This small discrepancy can be due to the relatively large 
fraction of X-ray sources at z>3, in contrast with the 
strong drop at that redshift for the X-ray undetected 
sample. We also found that X-ray detected sources are 
in general brighter than the X-ray undetected ones in 
both the 24 /mi and R bands, as can be seen in Fig.[4j A 
K-S comparing the i?-band magnitudes for the X-ray de- 
tected and undetected sources found that the hypothesis 
that they are drawn from the same parent distribution 
is rejected with >99% confidence. This difference in op- 
tical magnitude explains why the X-ray detected sample 
extends to higher redshifts. Furthermore, as mentioned 
before, the redshift completeness for the IR-red excess 
sources not detected in X-rays is 49%, while for the X- 
ray detected sources is 83%, thus supporting the idea 
that the optical brightness explains the difference in red- 
shift distribution. 

Compared to the general distribution of optically- 
detected sources in the ECDF-S, the IR-red excess 
sources are found at significantly higher redshifts. 



This may be an effect of the selection, namely the 
/(24/im)//(i?)>1000 requirement. Because for obscured 
AGN the rest-frame optical/UV light is dominated by 
the host galaxy (e.g.. ITreister et al.ll2005l and references 



therein), the / (24/tm) / f (R) is essentially a measurement 
of the AGN-to-host galaxy luminosity ratio. Hence, the 
fact that the 24 /(in-selected sources are found mainly at 
z>\ is a consequence of the available volume for high- 
luminosity sources. Furthermore, these sources do not 
have local analogs since even the most luminous nearby 
star-for ming galaxies ca nnot reach such high f^A^ml fn 
ratios (|Dev et aUfeOOSf K The high-redshift dropoff at 
z~3 can be explained by the faintness of the optical coun- 
terparts at such redshifts, thus making the determination 
of photometric redshifts very hard. 

5.2. X-ray Stacking 

Thanks to its very low background, Chandra data can 
be efficiently stacked in order to study the average X- 
ray properties of a sample of individually-undetected 
sources, thus achieving an effective exposu re time of 
many megaseconds (e.g.. iBrandt et afll2001l ). In order 
to perform X-ray stacking we used both the web-based 
CSTACK code 16 developed by T Miyaji and our own 
scripts, which are based on CIAO (jFruscione et al.ll2006h 
version 4.1. In all cases we found consistent results using 
these two codes. Stacking was done using the final un- 
binned mosaics of ECDF-S Ch andra observations, which 
were produced as described bv lVirani et a l. (2006). The 
stacking was performed in two bands independently: soft 
(0.5-2 keV) and hard (2-8 keV). 

We started from the sample of 185 24 /im-selected 
sources in the ECDF-S not detected in X-rays. To avoid 
possible contamination and enable a more accurate back- 
ground measurement, we excluded 46 sources for which 
there was a detected X-ray source closer than 15". Af- 
ter coadding the remaining sources the effective exposure 
time was ~30 Msec or ^1 year. The stacked image in 
both the soft and hard bands is shown in Fig. [6l To- 
tal source counts were measured in a 3"-radius circle. 
The background properties were estimated in two differ- 
ent ways: by using an annulus with a 7" inner radius 
and 15" outer radius and by performing photometry in 
random 3" circles outside the central region. The re- 
sults from both methods are roughly consistent with each 
other. We measured a total of 68 and 89 background- 
subtracted counts in the soft and hard bands. The back- 
ground standard deviations are 14.1 and 32.8 counts re- 
spectively, which corresponds to a signal to noise of ^4.8 
in the soft band and ~2.6 in the hard band. The signifi- 
cance of these detections was established by performing 
500 independent Monte Carlo simulations, stacking 139 
random positions in the field, thus having the same effec- 
tive exposure time and noise properties as the stacking of 
the 24 /im-selected sources. By studying the distribution 
of central background-subtracted counts in these simu- 
lations we found that, as expected, they are consistent 
with a Gaussian distribution with standard deviations 
of 17 and 22 counts in the soft and hard bands respec- 
tively. This indicates that the signal detected for the 
24-/im-selected sources is significant at the ^4-er level 
in each band. The latter measurement is more robust 

16 http://cstack.ucsd.edu/ 
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Fig. 6. — X-ray stacking of 129 sources not detected individually 
by Chandra in the soft (0.5-2 keV; left panel) and hard (2-8 keV; 
right panel) bands. The equivalent exposure time is ~30 Msec. A 
significant detection, ~4-<t, is found in each band independently. 
The total background-subtracted counts are 68.1 in the soft band 
and 86.5 in the hard band, which correspond to average fluxes of 
~2 and ~8 xl0 — 17 erg cm - 2 s — 1 respectively. 



because it is based on 500 independent determinations, 
compared to the twenty 3" circles used in the method 
described previously, hence significantly reducing the ef- 
fects of the faint undetected X-ray sources in the distri- 
bution of background counts. 

In order to convert count rates into fluxes we 
used the Portable, Interactive Multi-Mission Simula- 
tor (PIMMS) tool 17 . Assuming the corresponding 
Chandra/ACIS response functions, an intrinsic power- 
law spectru m with r=1.9. G alactic absorption of 
8xl0 19 cm -2 (jKalberla et al.ll2005l ). and intrinsic absorp- 
tion of 10 24 cm~ 2 at z—2, we found conversion factors 
from counts/second to erg cm _2 s _1 of 6.5xl0 -12 in the 
soft band and 2.6xlO -11 for the hard band. Hence, 
the average observed X-ray fluxes for the 24 /im-selected 
sources are ~2.1xl0 -17 erg cm' 2 s _1 in the soft band 
and ~8xl0~ 17 erg cm _2 s _1 in the hard band. If instead 
a lower column density, 2xl0 23 cm~ 2 , is assumed, the 
fluxes in the soft and hard bands are reduced by 17% 
and 20% respectively. Similarly, if a z—l is assumed in- 
stead for a NH—W 24 cm^ 2 spectrum the conversion for 
the hard band is increased by 27%, while if z=3 is used 
the conversion is decreased by 13%. We can hence con- 
clude that the uncertainties due to the assumption of an 
average spectral shape for the counts-to-flux conversion 
generate uncertainties smaller than ~30%. 

The HR for the stacked X-ray signal is 0.13±0.06. For 
comparison, for the sources individually detected in X- 
rays we found an average HR of -0.23±0.01, significantly 
softer than the value for the X-ray-undetected sample. 
This higher HR suggests the presence of more obscured 
sources in the X-ray undetected sample. Converting di- 
rectly this HR into a Njj for a intrinsic power-law spec- 
trum with r=1.9 and the average redshift of our sample, 
z~2, we find that N H ~1.8xl0 23 cm" 2 . If instead the 
extreme values of z—l and z=3 are considered we obtain 
values of Nh—5x 10 22 and 4x 10 23 cm -2 respectively, i.e., 
a uncertainty of ~2x can be expected due to the redshift 
distribution of our sources. However, and more impor- 
tantly, we note that the observed HR could be signifi- 
cantly affected by the presence of X-ray emission from 
non-AGN galaxies, which have in general a softer spec- 

17 Web version availa ble at 

http://heasarc.gsfc.nasa.gov/Tools/w3pimms.html 



trum (e.g. jFab biano 1989j). In order to estimate the frac- 
tion of obscured A GN relative t o star - forming galaxies in 
a stacked sample, iFiore et al.l ([2009.) performed simula- 
tions combining these two types of sources to match the 
observed HR. These simulations take into account the 
redshift distribution of our sources. Converting our ob- 
served HR into the bands used in their work (0.3-1.5 and 
1.5-6 keV) we find that the fraction of heavily obscured 
AGN with A^>10 23 cm - 2 is -90%. 

To understand the effects of the sources without either 
a photometric or spectroscopic redshift determination on 
our results we stacked separately the 90 sources with a 
redshift measurement. A total of 25 and 31 background- 
subtracted counts were detected in the soft and hard 
band respectively. The HR derived from these detections 
is 0.11±0.16, i.e., in very good agreement with the value 
found for the whole sample. This indicates that while 
sources without a measured redshift can have a different 
redshift distribution, this does not affect our results sig- 
nificantly. Similarly, splitting the sample with redshift 
measurements at z—2, roughly the median of the distri- 
bution, we detected 18 (soft band) and 27 (hard band) 
background-subtracted counts for sources at z<2, while 
no significant detection was obtained for the z>2 sample. 
The HR for the z<2 sample is 0.2±0.18, slightly larger 
but consistent with the value obtained for the whole sam- 
ple. 

To further study the possible effects of the spread in 
redshift of our sample we also performed X-ray stacking 
in the approximate rest-frame, i.e., multiplying each ob- 
served photon energy by 1+z before stacking. (This is 
not strictly correct since one should correct the incident 
photon energy rather than the observed photon energy. 
However, for the Chandra ACIS detectors, as for most 
X-ray detectors with modest energy resolution, a delta 
function of incident photons is redistributed over a broad 
range of energies and the response matrix cannot be in- 
verted in a stable or unique way.) This pseudo-rest-frame 
stacking can be done only for the 90 sources with redshift 
measurements, which is roughly half the sample. Exclud- 
ing the sources closer than 15" to a individually-detected 
X-ray source further reduces our sample to 63 sources. 
We perform this X-ray stacking in four rest-frame bands: 
1.5-5, 5-10, 10-15 and 15-24 keV. In the first three 
bands we obtained a signal-to-noise ratio greater than 
~2, while there is no significant detection in the last 
band, most likely due to the rapid decrease in sensitivity 
of Chandra at high energies and the reduced number of 
sources at high redshift in our sample. These results are 
shown in Fig. In order to compare with sources with 
a known spectral shape, we also stacked the X-ray emis- 
sion from the 15 IR-red excess sources individually de- 
tected in X-rays with a measured redshift. As described 
in section |4~T1 these sources are obscured AGN with an 
average A^~10 23 cm -2 . We also compare with the rest- 
frame X-ray spectrum of a luminous unobscured AGN 
found in the s ame fi eld, XID 312, which according to 
iTreister et all (|2009bD has a redshift of 1.887 and hence 
covers the same rest-frame energy range as our IR-red 
excess sample. As expected, the X-ray spectrum of the 
IR-red excess sources is significantly harder than those of 
the X-ray detected sources with similar optical/IR colors 
and a typical unobscured AGN. 

To quantify the amount of absorption from the rest- 
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Fig. 7. — Spectrum resulting from X-ray stacking performed 
in the approximate rest-frame. The black points show the 
background-subtracted counts measured from the combination of 
63 IR-red excess sources with measured redshifts that were individ- 
ually undetected in X-rays, while the red points show the combined 
signal from the 15 X-ray detected IR-red excess sources. The blue 
points show the spectrum of source XID-312, a luminous unob- 
scured AGN in the ECDF-S at z=1.887. The last two spectra were 
normalized to the measurement for the X-ray undetected sources 
at 10—15 keV, where the emission should be more isotropic. The 
green and brown data points show the expected spectral shape for 
heavily obscured AGN with column densities of JVj/ = 10 and 
>10 24 cm -2 (see text for details). The magenta data points show 
the combined spectrum of a population of 90% CT AGN and 10% 
pure star-forming galaxies. The spectrum of the X-ray undetected 
sources can be well-explained by either an obscured AGN with 
N H >W 23 cm- 2 or a mixed population of Compton-thick AGN and 
star-forming galaxies. 



frame X-ray stacked signal presented in Fig. we com- 
pute the X-ray spectra of a heavily-obscured AGN with 
Ay = 10 23 and >10 24 cm~ 2 . For the latter, we incorpo- 
rate the eff e cts of Compton scattering, as described by 
iMatt et alJ l| 19991 ). and an Fe Ka line at 6.4 keV with 
an eq uivalent width of 2 keV, typical of Compton-thick 
AGN (jGilli et al.lll999f ). Fig. shows that for the sam- 
ple of 15 IR-red excess sources individually detected in 
X-rays, an Njj of 10 23 cm~ 2 explains the observed ratio of 
5-10 to 10-15 keV flux; at lower energies, the contribu- 
tion from less-obscured sources in this sample can explain 
the observed excess relative to the model predictions. In 
the case of the IR-red excess sources individually unde- 
tected in X-rays, the spectrum in the 1.5-15 keV range 
can be well explained both by an obscured AGN with 
A#>10 23 cm~ 2 or by a mixed population of a majority 
of Compton-thick AGN and star-forming galaxies, with 
the later providing most of the flux at E<5 keV. For ex- 
ample, in Fig. [7] we show the resulting X-ray spectrum of 
combining a population of 90% CT AGN with 10% pure 
star-forming galaxies. This is fully consistent with our 
results derived from the X-ray signal stacked in the ob- 
served frame presented above. Hence, we conclude that 
the redshift spread does not affect our conclusions and 
therefore in this work we focus on the observed-frame 
measurements in order to increase our sample and thus 



the significance of our results. 

A possible dependence of the fraction of AGN on mid- 
IR flux (and thu s luminosity) for X-ra y sel ected sources 
was re ported by iTreister et al.l (|2006l ) and iBrand et al.1 
(2006). Hence, it is reasonable to expect a similar be- 
havior for our 24 ^m-selected sample. In order to test 
for this effect, we separated our sample in three bins 
based on 24 //m flux, with the number of sources chosen 
to guarantee a signal-to-noise ratio for the X-ray stacked 
signal greater than two. The brighter bin includes 34 
sources with 268< fe^m (A*Jy)<1292. The signal to 
noise of the stacked signal in the hard band is ~2.9. For 
these sources we meas ure a HR o f 0.37± 0.25, which us- 
ing the simulations of iFiore et al . (2009) corresponds to 
a AGN fraction of ~95%. At intermediate 24 /jm fluxes 
we stacked 33 sources with 155</24 /im (^(Jy)<266. The 
detection of the stacked signal in the hard band is 2.1. In 
this case we found a HR of -0.02±0.01, or an AGN frac- 
tion of ~80%. Similarly, for the faintest 72 sources, with 
53</ 24Mm (^Jy)<148 we found that HR=0.018±0.01, 
which also corresponds to an AGN fraction of ~80%, 
although in this case the hard-band detection is not sig- 
nificant (<7<2). In conclusion, although we find a weak 
trend of AGN fraction with 24 /im flux, it is not a sta- 
tistically significant result. 

5.3. Optical/ Mid- IR Spectral Fitting 

Because in obscured AGN the optical light is dom- 
inated by the host galaxy, it is feasible to study the 
stellar population in these sources. In order to do 
that, we fit the UBVRIzJK spectral energy distribu- 
tion of our IR-red excess sources to an ensemble of model 
star form ation histories based on a two-burst scenario 
(see e.g. iFerrerasfc Silkl 120021: ISchawinski et all 120071: 
iKavirai et al.ll2007t ) using the stellar models of lMarastonl 
( 20051) . We allow the young burst to make up between 
1% and 100% of the stellar mass and range freely in age 
up to the age of the old burst, which must be at least 
500 Myr old. We vary the decay timescale of the young 
burst from 10 Myr (an instantaneous burst) to 5 Gyr 
(effectively constant SFR) and fix the metallicity of both 
components at the solar value . We apply a scree n of 
dust using the extinction law of ICalzetti et alJ (|2000l ) for 
starburst galaxies and allow the color excess E(B-V) to 
vary from to 2.8. An AGN component was not added 
since for these sources the amount of obscuration of the 
nuclear region should be large enough to completely ab- 
sorb the UV and optical emission, in order to explain 
the red R-K color. Only those sources with detections 
in sufficient bands to have at least one degree of freedom 
(six bands) are fitted. For the majority of sources where 
we do not have a spectroscopic redshift, we assume the 
photometric redshift to be the true redshift. We find 
that not all IR-red excess sources result in an acceptable 
value of the reduced % 2 statistic: of those sources with 
X-ray counterparts, 12 out of the 15 that have detections 
in at least six bands and a redshift measurement, have a 
reduced x 2 < 3.84. 

The nature of the available data limits what we can 
infer about the star formation history of these sources. 
We are sampling the rest-frame UV and blue optical part 
of the spectrum and expect a large amount of extinction 
in the UV. We also only have broad-band photometry, 
so as a result, we cannot constrain whether the current 
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Fig. 8. — Spectral energy distributions for the 12 IR-red excess 
sources in the ECDF-S individually detected in X-rays with mea- 
sured spectroscopic or photometric redshift and detections in at 
least six bands. Circles with error bars show the observed photo- 
metric fluxes in the UBVRIzJHK bands. Dotted lines show the 
fit to the rest-frame optical light using galaxy templates. The me- 
dian stellar mass for this sample is 5x10 Mg, while the average 
extinction is E(B-V)=0.6. 



Fig. 9. — Spectral energy distributions for a random sample 
of 16 IR-red excess sources not detected in X-rays individually. 
Symbols are the same as in Fig. \E\ The mean stellar mass for 
the IR-red excess sources not detected in X-rays is 10 Mq, while 
the mean E(B-V) is 0.4. Both the stellar mass and extinction are 
slightly lower but consistent with the values found for the X-ray 
detected sources. In general we found evidence for a substantial 
young stellar population, indicating that the host galaxies of most 
of these sources are experimenting significant moderately-obscured 
star formations episodes. 



starburst that we do infer is an event in which most or 
all of the stellar mass is formed, or wether it is a less 
substantial burst with younger ages and higher amounts 
of extinction. In Fig. [5] we show the resulting UV-optical 
spectral fitting for the X-ray detected sources. The me- 
dian stellar mass for this sample is 4.6xlO n Mq, while 
the average extinction is E(B-V)=0.6. 

For our sample of IR-red excess sources not detected 
in X-rays, and requiring a redshift measurement and de- 
tection in at least six bands, we found acceptable fits 
for 60 out of 84 sources. The most likely explanation 
for the poor fits are inaccurate photometric redshifts. 
We marginalize over the nuisance parameters to investi- 
gate the stellar masses and levels of optical extinction in- 
ferred for our sample of IR-red excess sources. In general, 
they contain substantial stellar populations in an ongoing 
burst younger than 100 Myr and widely range in levels of 
optical extinction (0<E(B-V)<1) 18 . The best-fit stellar 
masses range between 10 9 -10 M© with a median stellar 
mass of ~10 11 M Q . We caution that the typical uncer- 
tainties in the individual stellar masses after marginal- 
ization are on the order of a factor of 10. We further 
investigate the effects of the uncertainty in the photomet- 
ric redshift derivation in our derived stellar masses and 
star- formation histories, by varying the redshift within 
its 1-sigma allowed range while doing spectral fitting. 
We found that the dispersion in the stellar masses is on 
the order of 0.1 dex, which is substantially smaller than 
the typical errors on individual masses. Furthermore, the 
typical properties of the derived star-formation histories 
are not affected. Hence, we conclude that for the sources 
with good spectral fits, the errors on the derived stellar 
masses and star-formation histories are dominated by er- 
rors on the fluxes, in particular in the near-IR bands, and 
not by the photometric redshift errors. 
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Fig. 10. — Left panel: Distribution of stellar masses for the IR-red 
excess sources individually detected in X-rays (hatched histogram) 
and not detected in X-rays (solid histogram). The median values 
and distribution for these two samples are very similar. Indeed, 
a K-S test performed to these distributions indicate that they are 
consistent with being drawn from the same parent distribution at 
the 90% confidence level. Right panel: Distribution of extinction 
in the host galaxy, as measured from template fitting to the rest- 
frame UV/optical light. Symbols are the same as in the left panel. 
Contrary to X-ray absorption, which measures mostly the AGN 
obscuration, no significant differences are found between these two 
samples. In summary, we did not find major differences between 
the host galaxies of the X-ray detected and undetected sources. 



In Fig. [5]we show examples of spectral fitting for IR-red 
excess sources not detected in X-rays. We plot the his- 
tograms of the distributions of best-fit stellar masses and 
optical extinctions in Figure fTUl - the hashed histogram in 
each panel represents the sources individually detected in 
X-rays, while the solid histogram represents those with- 
out detections. We perform a K-S to determine whether 
the E(B-V) and stellar mass distributions of the X-ray- 
detected and undetected IR-red excess sources are consis- 
tent with being drawn from the same parent population. 
We find that at the 90% confidence level, both distribu- 
tions are consistent, and therefore no significant differ- 
ences are found between the host galaxies of the X-ray 
detected and undetected sources. 



18 This range in E( B-V) corresponds to 0<Ay<4.05 assuming 
an R v = 4.05 that the Calzctti ct al. (2000) law uses. 
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6.1. Near and Mid-IR Colors 



As it was shown in Section 15.21 the stacked X-ray sig- 
nal for our 24 /zm-selected sample indicates that most of 
these sources are consistent with being heavily-obscured 
AGN. In addition to the X-ray spectral shape, the rest- 
frame near and mid-IR colors provide important clues 
about the obscuration levels in these sources and the 
importance of the AGN relative to the host galaxy. 
As expected from m ost dust re-emission models (e.g., 
INenkova et~aT1l2002ft . the AGN spectrum at ~3-10^m is 
significantly affected by the amount of material in the 
line-of-sight (due to the effects of self-absorption) , which 
should be related to the obscuration of the X-ray emis- 
sion. 

As can be seen in Fig. [TTJ the distributions of /^//s 
flux ratios for X-ray detected and undetected sources 
are significantly different. Indeed, performing a K-S 
test we found that the hypothesis that these distribu- 
tions were drawn from the same parent distribution is 
rejected at the >99.999% confidence level. The distri- 
bution of /24/ fa ratios indicate a much bluer mid-IR 
spectrum for the X-ray detected sources. The average 
values of $2a/ f$ are 8.7 for the X-ray sources and 24.9 
for the X-ray undetected sources. At the typical redshift 
of our sources, these bands trace emission at rest-frame 
~8 and ~3 m i crons. The dust re-emission models of 
INenkova et alj (|2008bt ) suggest that these observed val- 
ues could be explained by the same intrinsic torus param- 
eters and two different viewing angles. For example, for 
a torus model with N=10 (number of equatorial clouds), 
cr=30 (width of cloud's angular distribution), Y=30 (ra- 
tio of outer to inner radii), q=2 (exponent of cloud's ra- 
dial distribution assuming a power law) and 7V=40 (op- 
tical de pth of an individual c loud; parameter definition 
given bv lNenkova et al.ll2008al). wh ich are consistent with 
the models used bv lNenkova et a l. ( 2008b) to explain the 
mid-IR spectra of local Seyfert galaxies, the 8 to 3 mi- 
cron flux ratios are consistent with the observed values 
for viewing angles of 30° for the X-ray detected sam- 
ple and 90° (edge-on) for the X-ray undetected sources. 
These parameters are mentioned just as an example. By 
studying a total of 19,939 torus models spanning a wide 
range of parameters we found that in ~10% of them the 
resulting fi^l f% ratios are consistent with the observed 
values, taking into account the redshift distribution of 
our sources and assuming a viewing angle smaller than 
30° for the X-ray detected sources and greater than 70° 
for the undetected ones. 

As shown in the sample of ult ra-luminous IR galaxies 
presented bv lArmus et al.1 ([2007), in their Fig. 10, only a 
few sources reach values of the /24//s flux ratio as high as 
those observed in our sample of X-ray undetected IR-red 
excess sources. Furthermore, those sources are mostly 
classified as LINERs, while using cl assification schemes 
based on optical emission lines (e.g-. lKewlev et al.ll2006 l ) 
they are found in the "composite" region, i.e. they can 
be explain by a combination of an AGN and strong star- 
formation. Hence, we can conclude that the observed 
/24//8 ratios in our X-ray undetected sample strongly 
suggest the presence of a heavily obscured AGN in the 
majority of them. 

At the median redshift of our sample the Spitzer IRAC 
observations trace near-IR emission at ~1 to ^3 mi- 




F (24 /im)/F (8 /un) 



Fig. 11. — Upper panel: Ratio of 24 fim to R-band flux den- 
sity versus 24 /an to 8 fim flux ratio for the 24 ^tm-selected AGN 
candidates in the ECDF-S. Same symbols as in Fig. [T] are used. 
Lower panel: Distribution of 24 fim to 8 fim flux ratio for the X- 
ray detected (solid histogram, increased by lOx) and undetected 
(dashed histogram) sources. A clear difference between the two dis- 
tributions is seen, with X-ray detected sources having much bluer 
colors. 



crons. This spectral region is dominated by a combina- 
tion of stellar light, AGN continuum and hot dust emis- 
sion. There have been several studi es of the IRAC prop- 
erties of X-ray-selected AGN (e.g., iBarmby et al.l 120061 : 
ICardamone et af1l2008D. One of the main conclusions of 
the work of lCardamone et al.l (|2008f ) is that AGN show a 
large spread in near-IR colors an d thus the IRAC c olor- 
color AGN s electi on methods of lLacv et al.l (j2004f ) and 
IStern et all ((2005), which work very well to find lumi- 
nous optically-unobscured sources, are not particularly 
efficient in finding obscured/low-luminosity AGN. Fig. 1121 
shows two combinations of IRAC color-color diagrams for 
the X-ray sources and 24 /xm-selected AGN candidates 
in the ECDF-S. Interestingly, the X-ray detected sources 
with / 2 4 /.fR>1000 a nd R- K> 4.5 fall inside the A GN re- 
gions of lLacv et ail (|2004h and lStern et~afl (|2005l) . How- 
ever, the vast majority of the X-ray - undetected IR-red 
excess sources are found outs ide the IStern et al.l (2005) 
region and on the edge of the lLacv et al.1 (|2004l ) locus. 

As for the /^/js flux ratio, we also find in the IRAC 
colors significant differences between the X-ray detected 
and undetected samples. In the IRAC case, most of 
the differences are found at longer wavelengths. In fact, 
we found very small differences in the [3.6]-[4.5] color 
(average [3.6]-[4.5]=0.79 for X-ray sources and 0.69 for 
the X-ray undetected sample), while significant differ- 
ences can be seen in the [5.8]-[8] color (average of 0.99 
for X-ray sources versus 0.5 for the X-ray undetected 
sample). As at longer wavelengths, these differences can 
be well explained by the e ffects of self- abs o rption of the 
dust re-emission. Using the INenkova et all ((2008a) torus 
models, we found that for heavily absorbed sources the 
emission at wavelength shorter than ~5 /im is dominated 
by stellar light, as the AGN re-emission is severely ob- 



12 



Treister et al. 





Fig. 12. — Spitzer IRAC colors for the X-ray sources and ob- 
scured AGN candidates in the ECDF-S. Symbols are the same 
as in Fig, [T] Left panel: [3.6]- [4 .5] versus [5.8]-[8] co lors. Dashed 
lines show the Stern ct al. (2005) AGN region and the Brand ct al. 
(2009) threshold for PAH-dominated sources. The IR-red excess 
sources not detected in X-rays are mostly outside the AGN re- 
gion, mainly due to their relatively blue [5.8]-[8] color. Right panel: 
/s//4.5 versus /s.8//3.6 flux rations . The dashed region shows the 
AGN locus according to Lacy et al. (2004). As can be seen in both 
panels, most of the X-ray detected 24 /^m-selected sources are found 
in the AGN regions, suggesting that the rest-frame near-IR light 
is dominated by the nuclear emission. The fs/fi.5 flux ratio for 
the non- X-ray-detected sources is slightly lower than the one found 
in X-ray detected sources. Together with the difference found in 
Fig. [TT] suggests that these sources are more obscured versions of 
the X-ray detected AGN. 



scured. This explains why the 24 /im-selected sources 
not detected in X-rays fall in the same region of the 
[3.6]- [ 4.5] versus [5.81- [81 d iagram as inactive galaxies at 
z>2 (jBarmbv et al.ll2006l ). In the case of the X-ray- 
detected sources, the significant torus contribution in the 
observed-frame 8 (im band explains why these sources 
are displaced to the right in th is diagram, so that they 
fall inside the lStern etafl £2005) diagram. A very similar 
effect can be seen in the right panel of Fig. [T^l where self- 
absorption can explain the observed difference in /s//4.5 
flux ratio. 

The mid -IR colors of 24 /z m-excess sources were also 
studied bv lPope et all ((200cl in the GOODS-N field, to 
similar MIPS 24-/im depths. The first surprising re- 
sult is that they found that >90% of their sources with 
/24//ii>1000 also have R-K>4.5. In contra st, we found 
that t his fraction is 77% in our sample, while iFiore et al.l 
(2008) rep ort a fraction of ~ 50%. One explanation men- 
tioned by IPope et "ail {2008) is that the fainter 24 /im 
sources have in general blue R-K colors. However, in 
our sample we found that out of the 74 sources with 
/24//r>1000 and R-K<4.5 only 8 have / 2 4 Aim <100 //Jy. 
Another possibility is that the different K-band magni- 
tude li mit of each sample can explain these discrepancies. 
While IPope et all (|2008h do not quote their flux limit 
in the K band, they mention that a large fraction (75 
of 79, 95%) of the sources with /24/xm//i?>1000 are de- 
tected in the if-band. Similar ly, the if -band im aging in 
the GOODS-S region used bv lFiore et all (|2008| ) is deep 
enough to detected all their 24 /xm sources. Hence, it 
is unlikely that these differences are due to the if-band 
flux limit either. 

The main conclusion of IPope et "ail {2008?) is that a 
large majority, 80%, of their sources with /24//i?>1000 
(most of them also have R-K>A.5) are dominated by 
star-formation and not AGN activity. They base this 
conclusion in the low value of fs,/ fi.5 (their Figure 4), 
which is taken as evidence for a mid-IR emission domi- 
nated by a star-forming galaxy. In our sample, a large 
fraction of sources have /s//4.5<2 (their threshold for 
star-formation versus AGN), including most of the X- 



ray detected 24 /xm selected sources, in which the AGN 
nature derived from the X-ray luminosity and spectral 
shape is clear. We point out that these low /s//4.5 
values can easily be explained by the effects of self- 
absorption in the 8 (jm band, as can be seen using the 
iNenkova et al.l (l2008bf) to rus models, and as discussed 
above. IPope et al.l ( 20081 ) also use the X-ray spectral 
shape to support their conclusion of a low number of 
AGN in their 24 /zm-selected sample. However, as we 
study in m o re de tail in section 15.21 and was shown by 
IFiore et al.l (|2008l ). proper simulations are required to 
convert the observed X-ray spectral shape from a stacked 
sample that includes both AGN and star-forming galax- 
ies into an AGN fraction. 

6.2. X-ray to Mid-IR Ratio and AGN Fraction 

Both X-ray and mid-IR are tracers of AGN activity, 
either direct emission in the former or re-radiation in the 
latter, and of recent star formation, although the rela- 
tive X-ray emission is significantly lower in the latter. 
Hence, the X-ray to mid-IR ratio is particularly impor- 
tant in s eparating emission from AGN and star forma - 
tion (e.g.. lHorst et al.ll2006tlRamos Almeida et al.ll2007h . 
In Fig. [13] we show the hard (observed-frame 2-8 keV) 
X-ray versus mid-IR (at rest- frame 12 fim) luminosity 
ratio as a function of rest-frame 12 /j,m luminosity for 
the X-ray sources in the ECDF-S. The rest-frame 12 ptm 
fl uxes were obtained fro m the quasar template spectrum 
of iRichards et all (2006), normalizing it to the observed 
MIPS 24 /j,m flux, which at the typical redshifts of our 
sources is very similar to rest-frame 12 /im. This partic- 
ular wavelength was c hosen since it is near ly unaffected 
by obscuration (e.g., iTreister et al.l l2009bl ) and it was 
also used b y previous studies of the local AGN popu- 
lation (e.g.. iHorst et al.ll2006l ; iGandhi et al.ll2009D . thus 
allowing for a direct comparison with our high-redshift 
sample. Typical values of the conversion from observed- 
frame 24 /j,m to rest-frame 12/im are ~0. 91-1. 26, with 
an average of 1.0. 

The average X-ray to mid-IR ratio for the X-ray 
sources is 0.65, while the median is 0.4. As can be seen in 
Fig. H3] this ratio for the X-ray detected AGN is spread 
over one order of magnitude. One advantage of this ratio 
is that, at least for Compton-thin absorption levels, it is 
independent of obscuration. In fact, insignificant differ- 
ences were found for the average X-ray to mid-IR ratio 
values for obscured and unobscured sources in the X-ray 
detected sample. However, for more obscured sources, 
this ratio is expected to decrease, due to the effects of 
absorption in the hard X-ray bands. For example, for 
the 24 /j,m-selected sources individually detected in X- 
rays, which typically have A^~10 23 cm~ 2 , the average 
X-ray to mid-IR ratio is 0.43. 

In order to model the effects of obscuration in the hard 
X-ray to mid-IR ratio for CT sources, both photoelec- 
tric absorption and Compton scattering need to be taken 
into account. We have done th is using the template X- 
ray spectra of lMatt et al.l {1999) , computed using Monte 
Carlo simulations. For a heavily obscured source with 
A^//>3xl0 24 cm~ 2 the observed ratio decreases by 1-2 
orders of magnitude. In the case of the IR-red excess 
sources which are not individually detected in X-rays we 
computed X-ray fluxes from the measured stacked sig- 
nal. For each source, the average X-ray flux measured in 
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Fig. 13. — Ratio of hard X-ray to mid-IR (at rest-frame 12 /im) 
versus rest-frame 12 (im luminosity. Red circles show the X-ray 
detected sources with HR>-0.3, while the blue circles mark the 
sources with HR<-0.3. Black triangles with error bars show the 
average ratio for the X-ray-detected sources in bins of 30 sources. 
Circles enclosed by squares identify the X-ray detected IR-red ex- 
cess sources. The dashed line at a ratio of ~0.6 shows the average 
value of Lx/L\2 M m for all the X-ray sources. The dotted line 
shows the relation between i ntrinsic X-ray and 12 /im luminos- 
ity for local Seyfert galaxies ([Gandhi et al.l [20091) . Green circles 
show the location of the most luminous rest-frame 12 /im X-ray 
undetected sources outside the selection region. X-ray undetected 
IR-red excess sources are shown by olive circles and upper limits. 
X-ray luminosities were derived from the stacking signal. Dashed 
lines marked "lxlO 24 ", "3xl0 24 " and "lxlO 25 " show the effects 
of X-ray obscuration in the observed Lx/^l2/im ratio. The thick 
solid, dot-dashed and dotted lines show the expected Lx / Li2fj,m 
for star-forming galaxies, considering different recipes to convert 
star formation rates into X-ray and IR luminosities (see text for 
details). 



24 fim flux bins as described in section [5721 was assigned. 
For the fainter 24 /im sources, for which no X-ray stacked 
signal was measured, the 1-a background fluctuations are 
used as upper limits. Very similar and consistent results 
are obtained if the X-ray stacking is done binning sources 
based on 24 fim luminosity instead. 

As can be seen in Fig. [T3J the X-ray to mid-IR ra- 
tio for the IR-red excess sources is about two orders of 
magnitude smaller than the average value for the X-ray- 
detected sample and for most sources falls in the range 
expected for obscuration of ^5xl0 24 to 10 25 cm -2 . A 
small fraction of the IR-red excess sample, 20%, have X- 
ray to mid-IR ratios lower than 7x 10 -3 . It would be very 
hard to explain such low ratios using obscuration, and 
hence it is more likely that these sources are not AGN 
but star-forming galaxies. On the other hand, the ob- 
served ratio for the IR-red excess sample is significantly 
larger than for the sources outside our selection region, 
i-e. /24//ii<1000 or R-K <4.5, even at similar rest-frame 
12 /jm luminosities. The observed X-ray to mid-IR ra- 
tio for IR-red excess sources at X^m^lO 45 erg s _1 is 
roughly 10 times larger than the values expected from the 
X-ray versus m id-IR luminos i ty for star-forming galax- 
ies reported by iDonlev et al.l (|2008f ) . Similar results are 
found if the relation between star formation rate and X- 



ray luminosity propose d bvlGilfanov et al.l (12004D . which 
is consistent with the iRanalli et al.l (|2003h correlation, 
is used together with the relationsh ip between IR lu - 
minosity and star formation rate of iKennicuttl (|1998l) . 
In order to transform rest-frame 12 firm into total IR 
luminosity we u sed th e observe d correlations g i ven b y 
iCharv fc Elbal (pOOll) and by iTakeuchi et all ((20051) . 
The obtained rest-frame 2-10 keV luminosities are con- 
verted into observed-frame 2-8 keV assuming <z>=2 
and the typ ical X-ray spectrum of high-mass X-ray bi- 
naries (e.g., iLutovinov et al.ir2005h . namely r=1.0 and 
a high energy cutoff, E c , of 20 keV. As can be seen in 
Fig. 1131 the X-ray to mid-IR ratio for the IR-red excess 
sources is significantly higher than the expected value for 
pure star-forming galaxies, thus confirming the A GN na- 
ture for the vast majority of our sample. 

To investigate a possible dependence on redshift, in 
Fig. [14] we show the X-ray to mid-IR ratio as a func- 
tion of redshift. As can be seen in this figure, the IR- 
red excess sources are systematically above the sample 
with similar IR luminosities and bluer /24 ^ m //i? colors. 
At the same time, the IR-red excess sources are located 
in the region expected for heavily obscured AGN with 
A#>3xl0 24 cm -2 at all redshifts. This confirm that 
our conclusions are not affected by the redshift spread 
in our sample. It is also interesting that the fraction of 
star-forming galaxies in the IR-red excess sample inferred 
from the X-ray to mid-IR ratio is consistent with the 
value found from t he comparison of th e observed HR with 
the simulations of iFiore et all (|2009l ). While some trend 
is apparent, the evidence for a possible dependence of the 
fraction of AGN-dominated sources on mid-IR luminos- 
ity is not conclusive in our sample. Quantitatively, the 
fraction of star-forming galaxies (measured using a X-ray 
to mid-IR ratio of 7xl0 -3 as the dividing point between 
AGN and star-forming galaxies) in the 10 43 < L^^m 
(erg s^KlO 44 luminosity range is 4/17=23.5±20%, con- 
sidering only statistical uncertainties, while for sources 
with Li2/ im >10 44 erg s" 1 is 6/64=9.4±6%. Furthermore, 
for sources with Li 2Aim <10 43 erg s' 1 is 8/9=89^%. 
The number of sources is small for lower luminosity 
sources (Li2 A i m <10 43 erg s _1 ), mainly due to a combina- 
tion of the MIPS flux limit and our selection method. 
Using the same threshold in X-ray to mid-IR ratio, 
the fraction of star-forming galaxies for sources with 
£i2/jm>10 44 erg s _1 and outside our obscured AGN selec- 
tion region is 80±7%, thus showing the clear differences 
in X-ray properties for our 24 /im selected sample, even 
at the same mid-IR luminosities. 

6.3. Space density of CT AGN 

Given the strong evidence presented here for the AGN 
nature of the majority of the IR-red excess sources, and 
the very high obscuration levels that these sources show, 
we can use this sample to study the properties of the 
CT AGN population at z>2. While CT AGN can be 
abundant at high redshift, z>2, they only represent 
a small contribution to the extragalacti c X-ray back- 
ground , roughly 1-2%, as was shown by iTreister et al.l 
(2009aj). Because of the current uncertainties in the mea- 
surements of this background radiation and degeneracies 
in the assumed models, this integral constraint cannot 
be used to infer the number of heavily obscured sources 
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Fig. 14. — Ratio of hard X-ray to mid-I R r atio as a function of 
redshift. Symbols are the same as in Fig. 1131 The dependence of 
this ratio on redshift for a constant obscuration value is shown by 
the dashed lines. As can be seen here, the IR-red excess sources 
(olive circles) are mostly located in the region expected for AGN 
with obscurations >5xl0 24 cm -2 at all redshifts and are systemat- 
ically above the sample of sources with similar IR luminosities but 
bluer /24/jm/f.R colors. Based on this figure we can conclude that 
the redshift spread of our sample does not affect our conclusions. 



Fig. 15. — Three estimates (thick lines) of the hard X-ray lu- 
minosity function of C ompton thick AGN at 2=0, as measured by 
ITreister efaLl fl2009aT ) (upper left panel), Delia Ceca et al. (2008a; 
lower left panel), and Yencho et al. (2009; right panel). Sym- 
bols with error bars show the measurements of the space density 
of Compton thick AGN as a function of hard X-ray luminosity, 
de-evolved to 2=0 following the corresponding evolution for each 
work (see text for details). Filled triangles: t his paper, fi ll ed cir - 
cles: [Treister et al. (2009a), filled •pen tagons: IFiore eTaU ( 2009), 
open square: Daddi et al. (2007), star. Alexander et al. (2008). 



at high redshift. While ITreister et"afl (|2009aD discuss 
a possible stronger evolution for heavi ly-obscured AGN 
relative to the less-obscured sources, fDelTa Ceca et al.l 
(2008a) concluded that the luminosity functions of X- 
ray-selected CT and less-obscured AGN were consistent 
with each other. These potentially contradictory results 
are studied in more detail here. 

We measured the space density of CT AGN using our 
sample of IR-red excess sources. Because of the potential 
incompleteness of our photometric redshift determina- 
tions we restrict our study to z<3. Furthermore, taking 
into account the flux limit of the MIPS 24 /im data in the 
ECDF-S, the space density was measured only at z<l.l 
for sources with L24 Mm <10 44 erg s . These choices min- 
imize the effects of incompleteness in our sample. How- 
ever, other sources of incompleteness, like the K-band 
magnitude limit, required for our color selection could 
be also important. X-ray luminosities were estimated 
from the stacked X-ray signal, as described in Section 

EJ 

In Fig. [15] we present our measurements of the space 
density as a function of luminosity (triangles) together 
with the expectations from three different z =Q CT 
AGN luminos i ty fun ction s: iDella Ceca et al] (|2008af l. 
ITreister et all (|2009al) and lYencho et all 1)2009^ . For the 
latter, the X-ray luminosity function for all AGN was 
converted into a CT AGN luminosity function assum- 
ing the luminosit y dependence of the obscured AGN 
fraction given by ITreister et al.l (|2009af ) , normalized to 
match the INTEGRAL measurements at low luminosi- 
ties. Except for th e lack of a decline at low luminosi- 
ties in the work of IDella Ceca et ail (|200 8a). the three 
curves are consistent with each other. In order to com- 



pare with these z—0 expectations, the ECDF-S measure- 
ments, together with other values obtained from the lit- 
erature were "de-evolved" using the corresponding evo- 
lution form for each luminosity function. The differences 
between observations and expectations must be due to 
the assumed evolution in the latter. While in all these 
works a luminosity-dependent density evolution was as- 
sumed, different par ameter values were obtain ed in each 
case. For example, IDella Ceca et al.l (l2008al) assumed 
a stron g evolution for all AGN, while ITreister et al] 
(|2009al) used the values estimated bv lUeda et al l (|2003f) . 
which give a much shallow er evolution. As c a n be se en in 
this figure, w h ile bo th the IDella Ceca et al.l (|2008ah and 
I Yencho et al.l (|2009l ) provide an acceptable description 
of the observation al data, the values obtained using the 
lUeda et al.l (|2003f ) evolution fall systematically above the 
z=Q curve. Hence, these results indicate that a similar 
shape of the luminosity function is found for CT AGN 
and less -obscured sources , in agr eement with the conclu- 
sions of IDella Ceca et aLl (|2008af ). 

In order to study in more detail the evolution of the 
CT AGN space density, in Fig. [TB] we present our mea- 
surements of the CT AGN space density as a function of 
redshift. A reasonable agreement, in particular at z<l, 
is found between both our observed values and others 
obtained from the literature and the lumi n osity function 
and evolution assumed bv lYencho et al.l (|2009h . How- 
ever, at higher redshifts and higher luminosities clear 
di screpancies are found. This was already pointed out 
bv lTreister et al] (|2009al ). who concluded that this differ- 
ence of a factor of 2-3 could most likely be due to either 
incompleteness in the Swift/BAT and INTEGRAL CT 
AGN samples at z=0 used to fix the luminosity func- 
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Fig. 16. — Space density of Compton thick AGN as a function 
of redshift. Filled triangles show the measurements described in 
this p aper. Squares: space d ensity from the work o f ITozzi et al.l 
(2006). Star. Measurement by Alexander ct al. (2008). Pentagons: 
Values reported by IFiore et al.l (120091') . Solid lines show the ex- 
pected space d ensity of Compton t hick AGN from the luminos- 
ity function of lYencho et al.l (120091), wi t h the overall normaliza- 
tion fixed to the results of Trcister ct al. (2009a), while the dashed 
lines show the e x pectati ons based on the luminosity function of 
IDella C cca ct al. (20083)- Red symbols show measurements and 
expectations for Lx >10 43 erg s — 1 sources, while the blue symbols 
are for Lx>10 44 erg s . While for the lower luminosity sources 
a good agreement is found between observations and expectations, 
higher luminosity sources lie well above the luminosity function. 
Furthermore, the new measurements seem to indicate as strong in- 
crease in the number of high-luminosity heavily-obscured sources 
at 2>2. 



tion normalization (because reflection-dominated AGN 
are missed) or to contamination by other types of sources 
in the observed values at high redshifts. However, after 
adding the measurements obtained using the IR-red ex- 
cess sources in the ECDF-S it appears not only that the 
systematic difference is still present but perhaps more 
importantly that there is a strong increase in the num- 
ber of CT AGN from z~1.7 to 2.4. 

As can be seen in Fig. [T3J most of the sources with 
only upper limits to their X-ray to mid-IR emission ra- 
tios have mid-IR luminosities lower than 10 44 erg s _1 . 
Hence, excluding these sources from our measurements 
of the space density of CT AGN as a function of red- 
shift presented in Fig. [16] does not change these values 
significantly. For sources with Lx>10 43 erg s _1 if those 
with only upper limits are excluded the space density is 
reduced by ^30%. In Fig. [T5] where we presented the 
space density of CT sources as a function of hard X-ray 
luminosity, the effects of excluding sources with upper 
limits are more relevant for the less-luminous sources. 
In the 43<log(Lx [erg s _1 ])<43.5, the space density is 
reduced by ^70% if upper limit sources are excluded. 
Similarly, in the 43.5<log(ix [erg s~ 1 ])<44, the space 
density should be reduced by ^60% by excluding the up- 
per limits from the sample, while at higher luminosities 
the effect is less than 10% and hence negligible. 

The strong evolution in the number of CT AGN at 



high redshift presented in Fig. [16] was not indicated by 
any existing luminosity function. It is unlikely that 
this evolution is due to selection effects, as results from 
different fully independent surveys and selection tech- 
niques are c ombined in Fig. I16[ namely X-ray selected 
sourc es from ITozzi et al.l ()2006| ). 24 /zm-selected sources 
from IFiore et al.l (|2009n and our work, and the sam- 
ple of CT AGN found using mid-IR spectroscopy by 
lAlexander et al.l ((2008). If confirmed, this result indi- 
cates a very large number of heavily obscured AGN at 
z~2.5 followed by a strong decline, by a factor of ~3, 
at z=1.5. This result can be inte rpreted in the contex t 
of the galaxy evolution models of iHopkins et al.l (2008), 
where quasar activity is driven by galaxy mergers and 
the supermassive black hole is initially completely sur- 
rounded by dust, before radiation pressure removes it 
and a "classical" unobscured quasar is visible. The con- 
sequences of this observed rapid evolution are beyond the 
scope of this work and will be studied in further detail 
in an upcoming paper. 

7. CONCLUSIONS 

We presented a study of heavily-obscured (CT) AGN 
candidates selected from their optical, near and mid- 
IR colors. These sources were selected by requiring 
/24Atm//ii>1000 and R-K>4.5. Using this selection 
method, we found a total of 211 infrared-red excess 
sources in the ECDF-S. Of these, 18 were individually- 
detected in the Chandra 250-ksec observations of this 
field. These sources are moderately-obscured, with 
7Vff=10 22 -10 23 cm~ 2 , according to their X-ray spectra. 
Two of these sources have Nh>W~ 24 cm -2 , and there- 
fore are classified as CT AGN, including one at z=4.65. 
The average redshift for the X-ray-detected sample is 
<z>^2.37 (median=1.85), in agreement with previous 
studies. By performing spectral fitting to the observed 
optical to near-IR fluxes of the X-ray-detected sources we 
found that the host galaxies have typical stellar masses 
of lO n - 7 M and are subject to moderate extinction, with 
<E(B-V)>=0.6. 

For the X-ray undetected sources, we found a slightly 
lower average redshift of <2>~1.9, based on photomet- 
ric measurements only. We found that in general the X- 
ray-detected sources are brighter in both the 24 /zm and 
optical bands compared to the undetected ones. Taking 
advantage of the low Chandra background, we performed 
X-ray stacking, finding a significant (>3rx) detection in 
both the soft and hard bands. The average X-ray flux 
of these sources is ~8xl0 -17 erg cm _2 s _1 , indicating 
that in order to detected these sources individually, ex- 
posure times of ~10 Msec with Chandra would be re- 
quired. Taken at face value, the observed average HR of 
0.13 corresponds to a Nh of ^2xl0 23 cm -2 . However, 
this can also be interpreted as the effects of combining a 
population of 90% CT AGN and 10% star-forming galax- 
ies, which in general have softer X-ray spectra. From 
this analysis we found marginal evidence for a small de- 
pendence of the fraction of AGN relative to star-forming 
galaxies on 24 fj,m flux, going from ~95% for the bright- 
est sources to ^80% at the lowest flux bin. The spectral 
fitting analysis performed to these sources indicate that 
in general there is evidence for substantial young stellar 
populations, younger than 100 Myrs. This suggests that 
these sources are simultaneously experimenting signifi- 
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cant star-formation and heavily-obscured AGN activity. 
We did not found significant differences in the stellar 
masses and extinction values for the host galaxies of the 
X-ray detected and undetected IR-red excess sources. 

The X-ray undetected IR-red excess sources have red- 
der f2i/fs nux ratios than those detected in X-rays. The 
average jW/Zs— 25 is hard to explain assuming the ob- 
served spectra of ultra-luminous infrared galaxies, how- 
ever it can be well-explained by dust re-radiation mod- 
els in which the effects of self- absorption are important. 
Similarly, while the IRAC colors of most of these sources 
are outside the typical AGN region, this can be explained 
by the effects of absorption and/or a dominance of the 
near-IR emission by the host galaxy. In contrast, the 
X-ray detected sources have IRAC colors consistent with 
those of less-obscured AGN. The X-ray to mid-IR ra- 
tio for the X-ray-detected sources is similar to the av- 
erage value for the overall X-ray population. For the 
X-ray-undetected sources this ratio is ~2 orders of mag- 
nitude lower than for the sample individually detected 
in X-rays. This ratio is consistent with obscuration lev- 
els of ^5xl0 24 to 10 25 cm~ 2 , while it is significantly 
larger than the relative X-ray emission expected from 
star-formation activity alone. Using a constant threshold 
in X-ray to mid-IR flux ratio, we found that the fraction 
of AGN in our sample should be ^80%, consistent with 
the value found from the analysis of the stacked X-ray 
signal. We also found a similar dependence of the AGN 
fraction on 24 fxm luminosity ranging from ^90% at the 
bright end to only ~10% for the faintest sources. 

Finally, we studied the space density of sources implied 
by our sample of 24 /xm-selected heavily-obscured AGN 



candidates. While at lower redshifts and luminosities 
we found a good agreement with population synthesis 
models and extrapolations of existing X-ray luminosity 
functions to higher obscuration, we found significant ex- 
cesses for high- luminosity sources at high redshifts. Fur- 
thermore, we found a strong evolution in the number of 
sources with Lx>10 44 erg s _1 from 2=1.5 to 2.5. This 
can be interpreted as evidence for a relatively short-lived 
heavily-obscured phase before the strong radiation pres- 
sure removes the surrounding dust and turns the sources 
into an unobscured quasar. 
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